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1. Introduction 

Suppose you know the locations of post offices or cell phone satellites, and you want to know what 
regions they serve. Or maybe you know the locations of atoms in a crystal, and you want to know 
what a fundamental region looks like. There are lots of reasons you might want to make a tiling 
around a given discrete set of points. A natural way to do it is with "Voronoi tessellations" — so 
natural, in fact, that it has been rediscovered numerous times over the years. 

On the other hand, if you have a tiling, you might want to decorate each tile with a few points 
to create or destroy symmetry. Or you might look at all the vertices of the tiling-points where 
three or more tiles meet-to extract combinatorial information. If you have a tiling, there are many 
ways to extract a point set from it. 

So, you can get tilings from point sets and point sets from tilings: doesn't this give you a way to 
associate point sets to point sets or tilings to tilings? Once you have a map from a class of objects 
back to itself, you can take a dynamical systems viewpoint to analyze the situation. In this paper 
we are going to do exactly that, with a new dynamical system based on the vertices of Voronoi 
tessellations. 

For those uninitiated with the Voronoi tessellation, we begin with its definition and then give 
the definition of our dynamical system. From there, the remainder of §1 is spent exploring the 
evolution of simple point sets, using these simplified examples to develop both the intuition and 
vocabulary needed for more interesting cases. In §2 we give a new proof of a theorem, first proved 
in [■]] , quantifying the growth in size of point sets over repeated iteration. Following that we will 
point out some interesting corollaries and give some estimates on the growth rate. We devote §3 
to discussing what questions interest us most from the dynamical systems viewpoint. 

For now, let's turn to the definitions. 

1.1. Our dynamical system. We start with a discrete point set P C M^, which we call the gen- 
erating set, the members of which we refer to as the generators. To avoid technical complications 
we typically assume that P is finite. The Voronoi polygon of a point p £ P, denoted V{p), is 
given by 

^{p) = {x G M^l Hp — x|| < Hp' — x|| for all p G P} . 
Simply put, the Voronoi polygon of p contains every point in the plane that is closer to p than to 
any other member of P, or is equidistant between p and a nearby generator point. The Voronoi 
tessellation of P is given by 

T(P) = {V{p) \peP}. 

A point set and its Voronoi tessellation are given in Figure 1. 

Given distinct p, p' G P, the sets V{p) and V{p') are not necessarily disjoint; in fact, the boundary 
of each Voronoi polygon is shared with other Voronoi polygons. If V{p) D V{p') is a line, ray, or 
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(a) (6) 

Figure 1. A point set P, it's Voronoi tessellation, and ^{P). 

line segment, we call it a Voronoi edge and denote it ep^pi. If the intersection of three or more 
Voronoi tiles is a point, we call that point a Voronoi vertex. The set of all Voronoi edges and 
vertices are given by £{T{P)) and V(T(P)), respectively. 

It is useful to notice that e^^p' always lies on the perpendicular bisector of the line between p and 
p' . This gives a method for constructing Voronoi diagrams: for p € P, sketch each perpendicular 
bisector between p and and another member of P. Then the Voronoi polygon V{p) is the intersection 
of all half-planes created by the perpendicular bisectors. It is also useful to notice that a Voronoi 
vertex is equidistant from the generator points of the tiles it is in. For proofs of these properties 
and a wealth of other information about the Voronoi tessellations, [7] is an excellent source. 

Now, the set V(T(P)) constitutes a point set in its own right, and so one might naturally wonder 
what its Voronoi tessellation looks like. And we need not stop there~the Voronoi tessellation of 
V(T(P)) will have a vertex set, too, so how does its Voronoi tessellation behave? We have the 
makings of a dynamical system on the set V (M?) of all finite point sets in the plane. 

Definition 1.1. Let P (M^). We define the Voronoi iteration of P to be 

u{P) = V{T{P)). 

For n = 2,3, ... we define the n-th Voronoi iteration of P recursively: 

zy"(P) = V (T (u^'-^P))) = V (i^"-^(P)) . 

Figure 1 (b) depicts the Voronoi tessellation of the six points depicted in 1 (a), and so the lighter 
seven points — the vertices of the tessellation — constitute the Voronoi iteration of the original set. 
Let's begin looking at the simplest cases this dynamical system has to offer. 

1.2. The really trivial cases. We take the convention of using \P\ to denote the cardinaUty of 

P. If |P| = 1, then the single Voronoi polygon is ah of M^. Since is just one big tile, it has no 
vertices, and so = for all n > 1. Then up the ante: suppose P = \p\,P2\- The Voronoi 

tessellation then fractures the plane into two half-planes, split along the perpendicular bisector of 
P\P2, the line segment joining p\ to p^. The vertex set, however, is still empty, and so = 

for all n > 1. 



A DYNAMICAL SYSTEM USING THE VORONOI TESSELLATION 



3 



There are two possibilities when P has three points: either ah three are cohinear, or they he on 
a circle. The former yields = for all n > 1 and in the latter, ^{P) is a single point. No 

matter how large P is, the special cases of collinearity and cocircularity always work out like this: 

Proposition 1.2. All the points in P are collinear iff z^"(P) = for all n > 1. All the points in 
P are cocircular iff \i'{P)\ = 1 and u'^{P) = for n > 1. 

The first fact follows naturally from the observation that if p, g G P such that Cp^q S £ (T(P)), then 
Cp^q lies on the perpendicular bisector of the line segment pq. When all the points are collinear, 
the perpendicular bisectors are all parallel and thus do not intersect to produce new vertices. 
Conversely, if ^{P) = then no Voronoi edges intersect, which implies all Vorono'i edges are 
parallel, and hence P must be collinear. 

To prove the cocircularity result, we introduce the concept of an empty circle. This is a circle 
whose interior does not contain any generator points and whose boundary passes through three 
or more generator points. Such a circle gets its name since it is "empty" of generator points. 
Proposition 1.2 follows immediately from the next Proposition. 

Proposition 1.3. A point q gM."^ is a Voronoi vertex in V (T(P)) if and only if it is the center of 
an empty circle} 

In addition to proving Proposition 1.2, this proposition also lets us quickly find the vertex set of 
T(P). To do so, pick a noncollinear triple p, q, r G P, and look at the unique circle passing through 
them. If this circle is empty, then place a vertex at its center. Once you have checked every triple, 
every vertex will be accounted for — and you never had to sketch an edge! 

So, we know what happens when the point set is really small, or collinear, or cocircular. Let's 
move on to... 

1.3. A slightly less trivial case. Suppose |P| =4. To discriminate between configurations that 
aren't collinear or cocircular, we require a more sophisticated vocabulary. We say that a subset 
j4 C is convex if, given any distinct x, y G A, the line segment xy is contained in A. (For 
practice, try proving that Voronoi polygons are convex!) The convex hull of a set P is defined 
to be the smallest convex set containing P. We write CH{P) to mean the convex hull of P, and 
write dCH{P) to restrict attention to its boundary. 

Definition 1.4. For p G P, we say p is on the boundary of P and write p £ Bd(P) if p £ 
P n dCH{P). Otherwise, we say p is in the interior of P and write p £ Int(P). 

In general it is clear that 

(1.5) |Int(P)| + |Bd(P)| = |P|. 

We say that two distinct points p,p' G Bd(P) are neighbors on the boundary of P if 
pp' C dCH{P) and no p" £ P — {p, p'} lies on pp' . For noncollinear finite point sets, any point on 
the boundary has two distinct neighbors. 

It is not surprising that the Voronoi tiles of boundary points would have infinite edges. We 
denote by <S'ir(T(P)) and £i{T{P)) the finite and infinite edges of T(P), respectively. 

Proposition 1.6. Assume not all p £ P are collinear, and let Cp^ £ £ (T{P)). Then Cp^pi £ 
£i (T(P)) if and only if p and p' are neighbors on the boundary of P. Moreover, Cp^pi £ £p (T(P)) 
if and only if p and p' are not neighbors on the boundary of P? 

^See [7], pg. 61. 
^Follows from [ ], pg. 59. 
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This then imphes 

(1.7) |£:HT(P))| = |Bd(P)|. 

Returning to the case of |P| = 4, let us assume that the points in P are neither collinear nor 
cocircular. We may have |Bd(P)| = 3 or 4, and we are going to prove that Figure 2 represents the 
only possible types of Voronoi' iterations. 



(a) (6) 

Figure 2. Configurations satisfying (a) |Bd(P)| = 3 and (b) |Bd(P)| = 4. 

First let's assume |Bd(P)| = 3. Since |Int(P)| + |Bd(P)| = |P| = 4, say Int(P) = {p}. Proposition 
1.6 tells us that each of the edges of its Voronoi polygon V{p) must be finite. As \P — = 3, 
there can be at most three edges; since V{p) must be bounded, there must be exactly three edges. 
So V{p) is a triangle, and the three infinite edges promised by Proposition 1.6 extend from its 
vertices. Since |2^(P)| = 3, we conclude that |i^^(P)| = 1 and u^{P) = for all n > 3. 

Now let's assume that |Bd(P)| = 4 and that P is neither collinear nor cocircular. We will prove 
that |j^(P)| = 2. By Proposition 1.3 we know that |i^(P)| > 1. To prove that |j^(P)| < 2, we use a 
trick that will come in handy again, during the proof of our main theorem. 

For q G i^(P), let p{q) be the degree of the vertex q: the number of edges touching q. Then we 
have p{q) > 3 for each q £ i^(P). Summing the degrees of all the vertices therefore is greater than 
or equal to 3 • |i/(P)|. Counting edges instead we see that each infinite edge touches exactly one 
vertex but each finite edge touches exactly two vertices. Thus 

3-|K^')l< E P('Z) = l^:/(T(P))|+2.|f^(T(P))|. 

Since there are only two pairs of points in P that are not neighbors. Proposition 1.6 implies 
|^^f(T(P))| < 2. Thus 

3-|i/(P)| < 4 + 2 -2 = 8, 

and so |^'(P)| < 2, proving that in this case, |z^(P)| = 2. Moreover we can conclude that v'^{P) = 
for all n > 1. 
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1.4. The rest of the cases are all hard. The case |P| = 4 turns out to be the last case for 
which z^"(P) is guaranteed to equal the null set for large enough n. Indeed, the case |P| = 5 has 
resisted our attempts to understand it! Figure 3 depicts configurations of larger cardinality where 
the point set does not iterate to a simpler configuration, which in turn enriches (and complicates) 
matters. We must suspend our case-by-case analysis at this point. 




2. How BIG IS u{P)7 

Having played with the dynamical system for a while, we, along with our colleagues Allison 
Edgren and Olivia Gillham, decided to explore the question of what happens to the size of 
as n goes to infinity. We had some success: a formula [3] (unpublished) that tells us exactly how 
many points will be in i^(-P)! We are going to share this result with you as soon as we describe one 
more essential piece of geometric information. 

2.1. Relating degeneracy & cocircularity. Compare the configuration given in Figure 3 (a) 
with those sketched in Figure 4. Each are configurations of five points, yet each has a distinctly 
different iteration. We already have the vocabulary necessary to distinguish the configuration 
sketched in Figure 4 (a) from the others: they differ with respect to the number of points on the 
boundary. The difference between the other two, however, is more subtle. 

A meticulous reader might notice that all of the vertices in Figure 3 (a) are of degree three while 
there is a vertex in 4 (b) with degree four. A Voronoi diagram with a vertex of degree greater than 
three is called degenerate; if all the vertices of the Voronoi diagram are degree three we call it 
non-degenerate. We can quantify this concept of degeneracy with the number Ic{P)'- 

Definition 2.1. Let {Ci, . . . , Ck} be the set of all empty circles of a point set P. The number of 
instances of cocircularity of P is given by 

k 

/,(P) = j;(|C,nP|-3). 

i=l 
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(a) (6) 
Figure 4. Two configurations of 5 points. 

In a non-degenerate point set, every vertex will be of degree three, and so every empty circle will 
intersect exactly three points in P. In such a case, we compute Ic{P) = 0. When a vertex q is 
of degree /c > 3, this implies that |Cq n P| = k. In this case, we say that q contributes {k — 3) 
instances of cocircularity. This argument proves the following: 

Proposition 2.2. Let p{q) denote the degree of a vertex q G t^{P)- Then 

q(^u{P) 

When given a degenerate configuration, one can obtain a non-degnerate point set by simply shifting 
any cocircular points by an arbitrarily small amount; further, when given a non-degenerate config- 
uration, it remains non-degenerate under sufficiently small perturbations. Hence, when studying 
or applying Voronoi tessellations it is customary to ignore degenerate configurations: the modified 
non-degenerate configuration is usually "close enough" to the original, and if the points are collected 
from real-world instruments we can never assume that our measurements are exact anyway. 

If only we could make such an assumption! In our system, computing Ic{P) is usually only 
possible if you already know where the vertices are and, as we shall see, incorporating degeneracy 
makes determining the long-term behavior much more difficult. But we cannot simply throw out 
these configurations; there exist situations, like that which is depicted in Figure 5, that ruin things 
for us. A non-degenerate point set may iterate to a degenerate one, and so assuming P is non- 
degenerate does not guarantee ^{P) is as well. 

2.2. The theorem on counting vertices. We are ready to state and prove the theorem discov- 
ered by the second author and his colleagues in [•]] (unpublished). The proof we present here, due 
only to the second author, is much simpler than the original. 

Theorem 2.3. If not all p ^ P are collinear, then 

\v{P)\=2-\P\-\m{P)\-h{P)-2. 
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(a) (6) 
Figure 5. Ic{P) = while Ic{i^{P)) = 1. 



Proof. We can use the trick from §1.3, summing the degrees of the vertices in ^{P). Since P is non- 
cohinear, every infinite edge intersects exactly one vertex and every finite edge intersects exactly 
two. Hence, recalling from §1.3 that \£j (T(P))| = |Bd(P)|, we get 

PiQ) = \£i (T(P))| +2.\£f (T(P))| = |Bd(P)| + 2 • \£f{T{P))\. 

Euler's formula for finite planar graphs^ tells us that V — E + F = 1. We can apply Euler's 
formula to T(P) by ignoring the infinite edges and infinite Voronoi' polygons. Doing this yields 
\u{P)\ - \£f (T(P))| + |Int(P)| = 1, and so \£f (T(P))| = \u[P)\ + |Int(P)| - 1. Thus 

P{q) = |Bd(P)| + 2 • (IKP)I + |Int(P)| - 1) . 

gei/(P) 

Recalling from §1.3 that \P\ = |Bd(P)| + |Int(P)|, we have 

Y p(g) = 2-|K^)|+2.|P|-|Bd(P)|-2. 

On the other hand, by Proposition 2.2, 

Y p{q)=3-W{P)\+UP). 

q&y{P) 

Thus 

2 • \v{P)\ + 2 • |P| - |Bd(P)| -2 = 3- \v{P)\ + Ic{P) 
and so, solving for |i/(P)|, we find 

IK^')! =2-|P|-|Bd(P)|-4(P)-2, 

as desired. □ 



^See [n], p. 5 
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From this and our preliminary analysis in §1, several interesting observations may be made: 

Corollary 2.4. 

(i) // \P\ < 5 then v'^{P) = for some n G N. 

(ii) // \P\ = 5 then either v"'{P) = for some n G N or = 5 for all n G N. 

(iii) If \v{P)\ > \P\ then \P\ > 5. 

(iv) // |Int(P)| = 2 and Ic{P) = 0, then \u{P)\ = \P\. 

Theorem 2.3 alone is not enough to predict the long-term behavior of our system, since computing 
requires that we can also compute |Bd(z^"~"'^(P))| and Ic{h''^~^{P)). We have been able to 
find bounds on the size of v'^{P), and we present one of interest next. 

2.3. Bounds on the size of v"'{P). With no assumptions on the geometry of v'^{P) we can make 
an upper bound on the size of which we conjecture to be sharp. 

Proposition 2.5. For all n> 0, we have |z^"(P)| < 2" (|P| - 5) + 5. 

Proof. By induction on n. For n = 1, plugging the inequalities |Bd(P)| > 3 and Ic{P) ^ into 
Theorem 2.3 immediately gives us what we want. Proceeding inductively, we get 

|i/"(P)| = 2 • |i/"-^(P)| - |Bd(i^"-i(P))| - h{v'^--^{P)) - 2 < 2 • \v'^-^{P)\ - 5 

< 2 • (2"-i • (|P| - 5) + 5) - 5 = 2" • (|P| - 5) + 5, 

which completes the proof. □ 

We have also been able to find an upper bound on the size of the boundary of ^{P)- 
Theorem 2.6. The number of points on the boundary of u{P) does not exceed |Bd(P)|. 
We shall need the following lemma: 

Lemma 2.7. The interior angle for any vertex in a Voronoi polygon is strictly less than vr. 

Proof. Let p £ P generate the Voronoi polygon V{p) G T(P). Any vertex of V{p) lays at the 
intersection of some two Voronoi edges e^y and e^y, generated by distinct p', p" G P. From 
the argument given in §1.3, we deduce that V{p) is convex, and hence immediately have that the 
interior angle made by these edges must be less than or equal to vr. Seeking a contradiction, assume 
that this angle is equal to vr. By the perpendicular bisector property from §1.2, we may obtain p' 
by reflecting p across Cp^p' and p" by reflecting p across . But then p' = p" , contradicting the 
assumption that p' and p" be distinct. □ 

Proof of Theorem 2. 6. If all q G i^(P) are coUinear, then the theorem is trivial. To prove the result 
for non-collinear vertex sets, we seek to show |Bd(z^(P))| < |i£'/(T(P))| by proving that every 
q G Bd (i^(P)) intersects at least one infinite Voronoi edge. This, with the equality \£j (T(P)) | = 
|Bd(P)| from Proposition 1.6, gives us what we want. 

Since q G Bd(z^(P)), we may find neighbors q', q" G Bd(z^(P)) of q on the convex hull of i^(P). 
Define H{q, q') to be the closed half-plane through q and q' such that J^(P) is contained in H(q, q'), 
and similarly define H{q, q"); these half-planes exist by the properties of the convex hull. We thus 
find that J^(P) is contained in the intersection H{q,q') n H{q,q"), which we denote H. 

Let V{p) G T(P) be such that q G V{p) and V{p) n (M^ - i?) / 0. Let e G £" (T(P)) be an edge 
of V{p) that touches q. If e is infinite, then we are done. If e terminates in a vertex, then by the 
convexity of H we have that e is contained in H. Let e' £ £ (T(P)) be the other edge of V{p) that 
touches q. By Lemma 2.7, the interior angle made at q must be strictly less than vr. Since V{p) 
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must have a non-trivial intersection with — H, we find that e' cannot be contained in H. Thus 
it cannot terminate in a vertex, and instead is infinite, as desired. □ 

Note that the size of the boundary does not steadily decrease to three in all cases. There seem to 
be configurations with boundaries whose sizes stay stable over many iterations. 

Theorems 2.3 and 2.6 combine to give a lower bound on the size of v^{P) in the generic situation 
where there are no instances of co circularity. 

Theorem 2.8. LetP^V{M?) andNeN. Suppose that for all n = 0,1,2, ...N -1, Idi^'^iP)) = 0. 
Then \u^{P)\ > 2^\P\ - (2^ - 1) (|Bd(P)| + 2). 

A point set P with |Int(P)| > 2 and with Voronoi iterations that never contain instances of 
cocircularity will have exponential growth on the order of 2". (If |Int(P)| < 2 the size does not 
increase). Since having instances of cocircularity is a rare occurrence, we conjecture that such a 
point set exists. We show a typical situation in Figure 6: a few iterates of a randomly-selected point 
set of size 9. As Theorem 2.3 predicts, we see the size going to 13, then 21 points; if there continues 
not to be any cocircularity the growth will escalate rapidly: the next four iterations contain 37, 69, 
133, and then 261 points, respectively. 

3. So what's next? Questions from the dynamical systems viewpoint 

In dynamical systems terminology, V (M?) is called the "phase space" of the system, and elements 
of V (M^) are "states" of the system that evolve over time according to the map u. The orbit or 
trajectory of a state P is defined to be the sequence {P,u{P),i^^{P),...} ([1] or [G] are good 
references). So far we have only talked about one property of an orbit: the growth rate of |z^"(P)| 
over time. But there are many other interesting questions we can ask about this dynamical system. 

3.1. What is the right topology? In order to use the machinery of modern dynamics, we ought 
to have a topology on V{Mp) making Voronoi iteration continuous, at least some of the time. We 
need this to study major dynamical features such as recurrence — how orbits return to open sets 
over time — or topological entropy — a measure of the tendency of the system to become disordered. 
Finding the right topology is of great importance, but it has proved to be an interesting problem 
in its own right. We will mention some of the subtlety here. 

To begin, notice that the scale of the diagrams in Figure 6 increases with each iteration. We 
want our topology to care about the shape of the Voronoi polygons rather than their size. Similar 
triangles produce similar Voronoi diagrams, so the topology ought to respect that. Let us be more 
precise. 

Definition 3.9. A similarity transformation t G Sim(2) is a map of the form 

t(x) = kUx + xq, 

where x,xo G M^, U is a 2 x 2 orthogonal matrix, and k G M"*". For P, Q ^ Vi^), we say P is 
similar to Q, written P ^ Q, if there exists t G Sim(2) such that t{P) = Q. 

A similarity transformation is composed of translations, rotations, reflections, or dilatations of 
Euclidean space. One can show that Sim(2) forms a group under the composition of functions 
and that the relation ~ is an equivalence relation. Similarity transformations play nice with our 
dynamical system: 

Theorem 3.10. For t G Sim(2), we have t{v{P)) = i>{t{P)). 
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Figure 6. A few iterations of a point set with no instances of cocircularity 

Proof. Let q G '^{P)- By Proposition 1.3, there is a unique empty circle Cq, centered at q. Since 
similarity transformations preserve circles, t maps Cq to an empty circle centered at t(q). Thus 
t{q) G i'{t{P)), and so t{i^{P)) C i^{t{P)). By similar logic, since t{P) is a point set and is a 
similarity transformation, we have 

t-\u{t{P)))Ci^{t~'{t{P)))=i.{P). 

Thus z/ {t{P)) C t (z^(P)), and so we have equality. □ 

So the right topology should identify similar point sets so that we aren' t really looking at V(M?) 
but rather the quotient of V{M?) under similarity. 

On the other hand, forgetting the question of similarity, we do have an intuitive idea of what it 
means for point sets P and Q in ^(M^) to be "close". In fact there are already metrics, for instance 
the Hausdorff metric, for this. The problem is that we need our metric to respect Voronoi iteration. 

Suppose we draw little e-balls around all the points of P, and discover that each ball contains 
exactly one point of Q. If so, it is very likely that I'^P) and z/((5) are close as well. (The obvious 
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exception is with point sets that have nontrivial cocircularity — ^jigghng the points a Httle bit destroys 
the cocircularity and thus produces extra vertices in the Voronoi iteration.) This idea for a metric 
is promising but has a serious flaw: it can't compare sets that aren't the same cardinality. Since 
there are sets of the same cardinality that iterate to sets of different cardinalities, it would be nice 
to have the ability to consider whether those iterates are close. 

Unfortunately, trying to measure distance between sets of different cardinalities opens Pandora's 
box. If several points of Q are inside the e-ball around a point of P, we can create all sorts of 
bizarre patterns in y{Q) that may be quite dissimilar to v{P). We find ourselves unsure how to 
resolve these issues. 

3.2. Are there periodic points? For /c G N, we say P is a period-A; point if P ~ v^{P). As a 
special case, if A; = 1 then P ~ i^(P) and we call P a fixed point. Despite our attempts so far, we 
have not been able to find any finite periodic point set P. 

Curiously, it's not difficult to find an infinite set P of period 1 or 2! In Figure 7 (a), the square 
lattice iterates to a shifted copy of itself and so has period 1; in Figure 7 (b), points evenly spaced 
on the diagonals y = itx iterate to points on the x- and y-axes, which then iterate back for a period 
of 2. In both examples one sees a high degree of cocircularity and we assume that this may be a 
key component in finding examples with other periods as well. 




(a) (6) 
Figure 7. Infinite configurations which are (a) period-1 and (b) period-2. 

Once a fixed (or periodic) point has been located, we may ask whether it is attracting, repelling, 
or hyperbolic. We look at "nearby" points (again the need for a well-defined notion of distance!) 
and look at what their orbits do. If all nearby points come closer and closer to P, it is attracting; 
if they all get pushed away we call it repelling. If there is a mixture it may be hyperbolic. There's 
some indication that the grid in Figure 7 (a) may be hyperbolic: we know it repels grids with 
"defects" like the one in Figure 8. However, if we shift an entire row of points up by a fixed 
amount, the orbit will be pulled back towards the orbit in 7 (a). 

3.3. Sensitive dependence on initial conditions. Even though we lack a metric on our phase 
space V (M^), there is evidence of sensitive dependence on initial conditions. For instance, if we 
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shift a point on the grid in Figure 7 (a), it introduces a defect in the Voronoi iteration, which we 
see in Figure 8. As the system evolves, that defect will grow to include ever larger regions of the 
plane until the original grid is no longer recognizable. 
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Figure 8. Introducing a defect into the square lattice of 7 (a). 

One can see this effect in finite configurations as well. Consider the set depicted in Figure 4 (b) . 
Since |z^(-P)| = 4 we know it is doomed to iterate to the empty set. However, if we were to move any 
one of the four cocircular points even the slightest amount, we would see the single vertex of degree 
four become two vertices of degree three. This five-point configuration has a chance of iterating to 
five-point sets indefinitely. 

3.4. Can we go backwards? A sensible question to ask is, what happens when we try to invert 
the system? Most discrete point sets in the plane are not the vertex set of a Voronoi diagram. So 
we must ask first when a preimage of a point set exists, i.e., a set P' such that z^(-P') = P- If 
there is a preimage, when is it unique? When it is not unique, how many preimages are possible, 
and is there a "best" one? 

Some literature exists about inverting Voronoi tessellations (see e.g. ['', ']), and we know rel- 
atively straightforward conditions that must be satisfied for a preimage to exist. What is more 
subtle in our situation is that we don't have a Voronoi diagram — we just have its vertices. Without 
the edge adjacencies the problem becomes richer. 

It may be that no finite point set exists that can be inverted indefinitely, which would be a 
shame. It would be nice if there was a subset of P(M^) that was invariant under Voronoi iteration 
in both forward and backwards time. We conjecture that there is a subset S (ZV{^^)oi five-point 
configurations with 2 in the interior and 3 on the boundary for which (1) the set is invariant under 
forward Voronoi iteration, and (2) there exist preimages of all orders within the set. If we could 
identify criteria for membership in 5, then we could begin to examine the orbits of S for evidence 
of chaos. 

3.5. Conclusion. We admit that we don't know much about the behavior of this dynamical system 
yet. We know exactly what happens for small point sets, but as soon as there are five or more points 
we encounter difficulties. We have a theorem that measures how the point sets grow or shrink 
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in size, but it requires geometric information that isn't always readily available. We have evidence 
that many point sets grow without bound, but have been unable to determine which conditions 
guarantee this. We've noticed evidence of sensitive dependence on initial conditions, but we don't 
have the machinery to measure the phenomenon. In the course of our study we've developed tools 
to help us with our insight on this finicky dynamical system, and we've shared some of that here. 
There is one thing we know for certain: plenty of problems remain. Some are simple enough to be 
studied by budding mathematicians, and some may be subtle enough to interest their advisers as 
well. 
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